Fermi polarons in two dimensions 
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We theoretically analyze inverse radiofrequency (rf) spectroscopy experiments in two- component 
Fermi gases. We consider a small number of impurity atoms interacting strongly with a bath of 
majority atoms. In two-dimensional geometries we find that the main features of the rf spectrum 
correspond to an attractive polaron and a metastable repulsive polaron. Our results suggest that the 
attractive polaron has been observed in a recent experiment [Phys. Rev. Lett. 106, 105301 (2011)]. 
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The behavior of a mobile impurity (polaron) inter- 
acting strongly with a bath of particles is one of the 
basic many-body problems studied in condensed mat- 
ter physics With the advent of ultracold atomic 
gases [5|, the Fermi polaron problem in which a single 
spin-^ atom interacts strongly with a Fermi sea of spin- 
t atoms, has become a subject of intensive research [G^. 
In three dimensions it was found that the polaron state 
splits into two branches, a low-energy state interacting 
attractively with the bath of fermions, and the repul- 
sive polaron, which is an excited, metastable state 0- 
[9]. In this way the polaron exemplifies a more general 
paradigm of a many-body system driven into a nonequi- 
librium state where a small number of high energy ex- 
citations interact strongly with the surrounding degrees 
of freedom 0, [nj. The polaron is the limitine; case of 
a Fermi gas with strong spin imbalance, and the repul- 
sive polaron provides insight into the question whether a 
quenched, repulsive Fermi gas may undergo a transition 
to a ferromagnetic state even though it is highly excited 
0-[2l, [l3, [HI . Similarly, the ground state of the polaron 
problem has important implications for the phase dia- 
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gram of a strongly interacting Fermi gas 

It is a key question how many-body properties are af- 
fected by reduced dimensionality, and the polaron is a 
case in point. The combination of optical lattices and 
Feshbach resonances [5] provides a unique setting to ex- 
perimentally study strongly interacting low dimensional 
systems using ultracold atoms [17, 18]. Recent advances 
in radiofrequency (rf) spectroscopy afford to measure en- 
ergy spectra [14] and give access to excited states as 
well as full spectral functions using momentum resolved 
rf [H, So far, only the ground state of the two- 

dimensional polaron problem has been investigated the- 
oretically [21-23] with the focus on a possible polaron to 
molecule transition. This is similar to the 3D situation 
where for strong interactions it becomes energetically fa- 
vorable for the impurity to form a molecular bound state 
(2^ . The structure of high energy excitations and the 
experimental polaron signatures in rf spectroscopy have 
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FIG. 1: (color online). Left panels: spectral function 
Ai{q^ E — /Ji) for impurity atoms interacting with a 2D Fermi 
sea. Red lines indicate the free particle dispersion and white 
(black) dashed lines mark the dispersion of the attractive (re- 
pulsive) polaron. Right panels: corresponding rf spectra 
illustrating how weight is shifted from the attractive polaron 
state (peak at negative frequencies) to the new repulsive po- 
laron state at positive frequencies. The two-body bound state 
energy is (a) sb/sf =0.1, (b) sb/sf — 1, (c) £b/sf = 5. 



remained open questions which we address in this work. 
We derive the spectral functions of both the molecule 
and the impurity atom (Fig. [T] left) and find that the 
impurity state splits into the attractive and the repul- 
sive branch. We compute rf spectra for homogeneous 2D 
systems (Fig. [1] right) as well as for the experimentally 



2 



relevant quasi-2D geometries (Fig. |4]). Finally, we argue 
that our calculation provides an alternative explanation 
of the recent experiment by Frohlich et al. [l7| in terms 
of the polaron picture. 

A quasi-2D geometry can be realized experimentally 
using an optical lattice in one direction with associated 
trapping frequency ujz- In this case, a confinement in- 
duced two-body bound state exists for an arbitrarily 
weak attractive interaction 25-23 with binding energy 
£b > 0. The spatial extent of the bound state is related 
to the 2D scattering length given by a^D = ^/ ^Jtosb > 0. 
In the weak coupling BCS regime of small 3D scattering 
length < [5] these dimers are large and weakly 
bound {sb ^ ^z)] in the BEC limit of small asB > 0, 
the weakly interacting molecules are too tightly bound to 
feel the confinement {sb ^ ^^/(^<^3d) ^ ^z)- Around 
the Feshbach resonance {a'^^ = 0) there is a strong cou- 
pling regime where the binding energy attains the uni- 
versal value Sb = 0.244: huOz j5l.l27|. 

At finite densities the majority atoms form a Fermi 
gas with Fermi energy sf and the two-body scattering 
is replaced by many-body scattering which gives rise 
to important qualitative differences, most notably the 
emergence of two polaron branches. Spectral weight is 
shifted from the attractive to the repulsive polaron in the 
non-perturbative regime where the interaction parameter 
1/ ln(£^/2£i?) diverge s 12 8| and the confinement induced 
resonance appears [5|, |29| . 

We consider a two-component 2D Fermi gas in the limit 
of extreme spin imbalance, described by the grand canon- 
ical Hamiltonian 



ka 



E 

kk'q 



with single-particle energies Ska = /2ma- for species a 
{h= 1), chemical potentials jia and system area A. Hav- 
ing in mind the experiment of Ref. [17] , we focus on the 
case of equal masses = = m. Generalizations to 
mass imbalanced situations are straightforward 21, 22| . 
In the low-energy limit the attractive 5-wave contact in- 
teraction g can act only between different species due to 
the Pauli principle. The majority atoms are not renor- 
malized by the presence of a single impurity with finite 
mass such that ji^ = Sp = k'^/2m at zero temperature. 
The chemical potential /i^ of the impurity atom is de- 
termined such that the impurity state ||) has vanishing 
macroscopic occupation. Furthermore, fi^ is negative due 
to the attractive interaction between t and | atoms. 

Dressed molecule. The two-body scattering of a spin- 
t atom and a spin-^ atom is described by the exact two- 
body T-matrix [26| 



ToiE) = 



47r/m 



\n{eB/E)^i7r 



(1) 



The pole of the T-matrix at E — —sb corresponds 
to the molecular bound state, and the associated vac- 
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FIG. 2: (color online). Molecular spectral function Amoi{Q, 
for different values of the two-body binding energy sb/sf- (a) 
0.1, (b) 0.5, (c) 1.0, (d) 2.5. The dashed lines mark the log 
continuum cjq, dash-dotted and solid the root continuum uj±. 



uum scattering amplitude for two particles with rela- 
tive momenta k and —k in the center-of-mass frame is 
f{k = |k|) = mTo{2ek) = iTr/Ml/Palj,) + in] [5]. 

In the presence of a Fermi sea of spin-^ atoms, the 
molecular state is dressed by fluctuations and described 
by the many-body T-matrix. This can^be calculated in 
the Nozieres-Schmitt-Rink approach [30|, as done in the 



2D case by Engelbrecht and Randeria [3l|, |32| . We gen- 
eralize these results to the case of spin imbalance and 
obtain 



T-\q,uj) = T-\uj + iO + /i^ + /i; - eq/2) 

(27r)2 + iO + /i^ + /i; - - eu+q ' 



(2) 



with the Fermi function npi^s). At zero temperature 
where fi^ = Sp and fi^ < 0, we obtain an analytical 
expression for the many-body T-matrix 



T(g, uj)=To(^^z±^^{z- eq)^ - AepSq 



(3) 



with z = uo^iO — Ep^fi^ and ± = sgnRe(z — £q). Due 
to the constant density of states in 2D the many-body T- 
matrix can be expressed as the two-body T-matrix with 
the argument shifted by Pauli blocking. The molecular 
spectral function AmoiiQ-,^) = — 2 ImT(gr, cj) is shown 
in Fig. [2] for several values of the interaction strength 
parametrized by the two-body binding energy Eb- One 
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observes a bound state peak at low energies and the 
particle-particle continuum at higher energies. 

The continuum of dissociated molecules arises mathe- 
matically from the branch cut of the square root ([3]) in 
the region uj-{q) < u < uj-^{q), iju± = £f(1 ± Q/kp)'^ — Mi 
(dash-dotted/solid lines), as well as from the branch 
cut of the logarithm ([T]) for > 0J-\-{q) and for ujq = 
eq/2 — £f — IJi^ < 00 < uj-{q) \i q> 2kF (dashed lines). 

The bound state pole of the many-body T-matrix has 
the dispersion relation (2l| 



gq/2(gq/2 - Sf) + gB(gF " ^b) 

eq/2 ^ SB 



(4) 



which changes qualitatively with the two-body binding 
energy sb (see Fig.[2j): ioi eb <2£f the bound state has 
minimum energy at a finite wave vector with positive 
effective mass vff jm = (2 — 2//cf<^2d)~"'^ [33]. 

Polaron and quasi-particle properties. The impurity 
atom is dressed with virtual molecule-hole excitations 



and becomes a quasi-particle with self-energy j32l-l3 



■T{k^q,ek-fi^^uj), (5) 



which leads to the same ground state energy as a vari- 
ational ansatz jsHj. Here we have used the fact that in 
the zero-temperature polaron problem the molecule has 
vanishing macroscopic occupation. Hence it has spectral 
weight only at positive frequencies (cf. Fig. [2j) where the 
Bose distribution vanishes. We perform the integral in 
([5]) numerically and obtain the spectral function of im- 
purity atoms 

Ai{q, co) = -2 lm[uj + iO + /i; - £q - S;(gr, uj)]-\ (6) 

The frequency and momentum dependence of the spec- 
tral function is shown in Fig. [T] (left panel) for three values 
of the interaction strength. In Fig.[3]we display the zero- 
momentum spectral function Ai{q = 0,^ — /i^) versus 
interaction parameter 77 = ln(A;i?a2D) = —\n{eB/2£F)/2. 
In both figures we set the reference energy to the free 
atom threshold by subtracting the chemical potential /i^. 

At weak binding Sb <^ £f (Fig- Ilk) the attractive po- 
laron is a well-defined quasi-particle at small momenta 
but for g > /cf it scatters off virtual molecules and ac- 
quires a large decay width. For intermediate binding 
(Fig. [TJd) a new repulsive polaron state appears at posi- 
tive energies. It is a metastable state with broad decay 
width, and it is shifted to higher energy due to the re- 
pulsive interaction with the Fermi sea of spin-^ atoms. 
The dispersion of the repulsive polaron has a minimum 
at finite momentum q ^ kF refiecting a similar feature in 
the molecular spectral function (Fig. [2t); for larger mo- 
menta it approaches the free particle dispersion. Finally, 
for strong binding (Fig. [It) both polaron branches are 
well separated and the repulsive polaron becomes an in- 
creasingly long-lived and stable quasi-particle. Between 



the attractive and the repulsive polaron branches appears 
the molecule- hole continuum (see also Fig. [3]). Its spec- 
tral weight is small in the case of a broad Feshbach res- 
onance studied here, but it is enhanced for narrow reso- 
nances by an admixture of closed channel molecules [H. 
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FIG. 3: (color online). Polaron spectral function Ai{q = 
0, ^ — /j^i) versus the interaction parameter 77. The dashed 
lines indicate the perturbation theory of Ref. 28]. Left inset: 
effective mass m*/m of the attractive and repulsive polaron 
as well as the molecule. Right inset: crossover of the quasi- 
particle weight Z from the repulsive to the attractive polaron. 



It is instructive to see how the quasi-particle prop- 
erties of the polaron change as the interaction param- 
eter 77 is varied. The right inset of Fig. [3] shows a con- 
tinuous crossover where the quasi-particle weight Z = 
l/[l — dujTi{q = 0, uj)] evaluated at the quasi-particle pole, 
shifts from the attractive to the repulsive polaron branch: 
for small binding (77 > 0), the attractive polaron is the 
dominant excitation and the weight is gradually trans- 
ferred toward the repulsive branch for increasing binding 
(77 < 0). This crossover is also refiected in the effective 
mass m*/m (Fig. [3] left inset). Our strong coupling cal- 
culation reproduces the perturbative results [28] for the 
attractive and repulsive polaron energies in the weak and 
strong binding limits (dashed lines in Fig. [3j). 

Radiofrequency spectroscopy. The spectral properties 
of the imbalanced Fermi gas can be accessed experimen- 
tally using rf spectroscopy. We assume that an rf pulse 
is used to drive atoms from an initial state \i) to an ini- 
tially empty final state |/). We choose the final state 
to be strongly interacting with a bath of a third species 
It) such that I/) is in fact the impurity state, |/) = W). 
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This inverse rf procedure interchanges the roles of \i) and 
I/) with respect to Ref. [14]; it has been proposed in 
Refs. sl, [qI and reahzed in the experiment by Frohhch et 
al. [ij 

Within hnear response theory, the rf transition rate is 
given by [36] 



/rf(^rf) = 



-21]?flmx''(- 



-CJrf - Mi + /i/) (7) 



where l^rf is the Rabi frequency, cjrf the detuning of 
the rf photon from the bare transition frequency and 
fii(f) the initial (final) state chemical potential. The re- 
tarded correlation function can be computed from the 
corresponding time-ordered correlation function —i6{t — 
0([V^}(r,t)V^,(r,t),V^|(r^tOV^/(r^m HE^. In gen- 
eral vertex corrections are crucial [38', '39^, but we find 
that they vanish in the case of negligible initial state 
interactions as appropriate for the experiment [l7|. At 
T = 0, we obtain 



^rf(^rf) = 



A^{q,uJri^eq - iii) . (8) 



The integral in equation (|8j) is calculated numerically and 
the resulting rf spectra are shown in Fig. [T] (right panel). 
The rf probes the final |^) state spectral function along 
the free-particle dispersion up to the initial state chemical 
potential fii. As in the experiment [l^, we assume a 
balanced initial state mixture with jii = ji^. We find 
a peak in the rf spectrum once the detuning cjrf reaches 
the final state chemical potential ji^ {ji^ is negative in the 
polaron problem). The transfer of spectral weight from 
the attractive to the repulsive polaron can be directly 
observed in Fig. [H 

Comparison to experiments. In order to relate our 
results to harmonically confined Fermi gases, we have 
to connect the strict 2D calculation to the quasi-2D ge- 
ometry relevant to experiments [5|, \^ . Well below the 
confinement energy hjjz where only the lowest transverse 
mode is occupied, this can be done by replacing £b with 
the exact quasi-2D two-body binding energy. Thus sb 
becomes a function of both the 3D scattering length 
and the confinement length 1^ = y/h/muOz ([l^]? cf. 
Eq. 82 in [El). 

Recently the quasi-2D geometry has been realized ex- 
perimentally with a Fermi gas of ^^K atoms [17]. Fol- 
lowing the inverse rf procedure described above, an ini- 
tially non-interacting balanced mixture is driven into a 
strongly interacting final state. As long as its occupation 
remains small the final state is a Fermi polaron, and our 
calculation predicts the experimental rf response. 

In Fig. m we show our trap averaged rf spectra versus 
magnetic field. We use the experimental parameters of 
Ref. [17] with o;^ = 27r x 80 kHz, uj± = 27r x 125 Hz, 
and express asD in terms of the magnetic field [HI, lioj . 
We incorporate the radial trapping in the 2D plane using 
the local density approximation; the local Fermi energy 



is Si? (r) = Si? — muj^r'^ /2 with peak Fermi energy Sf = 
9 kHz. Finally, we average over 30 pancakes in the z 
direction 17l . 
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FIG. 4: (color online). Trap averaged rf spectra of a quasi-2D 
Fermi gas: rf detuning versus magnetic field B. The experi- 
mental data points (blue and red diamonds) are taken from 
Ref. [17]. Also shown are the energy of the repulsive (dashed) 
and attractive (dash-dotted) polaron as well as the two-body 
binding energy (solid, white) in a homogeneous system. 

We observe that the lower branch of the experimen- 
tal spectra (circles) agrees well with the attractive po- 
laron picture (Fig. |4|) and our calculation provides an 
alternative interpretation to the two-body bound state 
(solid line) put forward in Ref. [17]. We note that also 
the measured frequency shift in 3D as shown in ^/Ti] fits 
the polaron picture [8]. Our results show a second rf 
peak at positive detunings corresponding to the repulsive 
polaron. The dashed line in Fig. 4 indicates its quasi- 
particle energy in the bulk (cf. Fig. [3]). As similar for 
the attractive polaron energy (dash-dotted line) the trap 
average leads to a significant shift of the rf peaks to lower 
energies. The experimental data (diamonds) in this mag- 
netic field range agrees qualitatively with our calculation. 
One possible reason for the remaining discrepancy is the 
large final state occupation in the experiment. 

In conclusion, we studied Fermi polarons in two dimen- 
sions which exhibit an attractive and repulsive branch 
and computed their rf spectra. Additional work is needed 
to understand discrepancies between theory and experi- 
ment for repulsive polarons. As an example, pump and 
probe experiments in the form of a sequence of two short 
pulses may shed further light on this issue. 
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Supplementary material 

Many-body T -matrix. The molecular spectral func- 
tion displayed in Fig. [2] changes qualitatively with the 
binding energy. For large binding sb > (Fig. [2]i), 
the bound state (j4]) has minimum energy = Sp — Sb — 
at g = and we find an effective mass m*/m |g=o = 
2/(1-2£f/^b) > 2. While it is well known that 
the T-matrix does not yield the correct binding en- 
ergy of the molecule in the BEC limit [22] the repul- 
sive polaron, which is the focus of our work, is repro- 
duced correctly [28|. For smaller binding sb < ^sp 
(Fig. [2]3-c), the effective mass at g = becomes nega- 
tive and a new minimun i ap pears at finite wave vector 
Qc = 2y^kFCL2B — 1/<^2D [21| and with positive effective 
mass m'^/m\q=q^ = (2 — 2/kFa2B)~^ - For decreasing 
binding energy sb the bound state eventually touches 
the continuum at momenta q±/kF = 1 ± ^J\ — 2eBl 
(dotted line in Fig.[2K). For q- < q < qj^ the bound state 
has negative residue and it ceases to exist. 

Polaron self- energy. We compute the many-body T- 
matrix in the ladder approximation which is represented 
diagrammatically in Fig. [SJi. The self-energy of the | 
atom is given by scattering an ^ hole off a molecule as 
depicted in Fig. [SJd. Explicitly, the self-energy reads 

_____ _ [nB{z)Gf{k-q, z-uj) lmT{k, z) 

- nriz) ImG^ik, z) T{k ^q,z^ u)] (9) 

where G^'^ refers to advanced/retarded t Green's func- 
tions. The first contribution comes from the pole and 
branch cuts of the T-matrix. In the polaron problem 
neither the ^ state nor the molecular state are macroscop- 
ically occupied at zero temperature, hence the molecular 
spectral function oc ImT has weight only at positive fre- 
quencies 2; > where nB{z) vanishes, cf. Fig.O The sec- 
ond contribution of ([9]) with the bare t spectral function 
ImG^(/c, z) = —i7rS{z — £fc + fi^) directly yields Eq. ([5]). 
The I self-energy yields the | Green's function via the 
Dyson equation (|6]) depicted diagrammatically in Fig.[5fc. 

Local density approximation. To incorporate the ra- 
dial trapping potential we use the local density approx- 
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FIG. 5: Diagrammatic representation of (a) the rf photon 
self-energy, (b) the self-energy of spin-^ atoms, (c) the Dyson 
equation for the dressed Green's function C^, and (d) the 
many-body T-matrix between states |t) and ||). 



distribution of the density of non-interacting fermions 
within a single layer. I^^\uj^r) is the rf response in 
Eq. dHJ computed for a homogeneous system with local 
Fermi energy £p\r) and interaction parameter r]^^\r) = 
— ln{eB/'^Sp\r))/2. In the experiment [17] 30 layers 
have been populated. To obtain the complete rf response 
of the trapped system we finally sum over 30 contribu- 
tions /trap where we assume that the peak density of each 

layer, e^^O), varies according to a Thomas- Fermi profile 
along the lattice direction. 



Shape and length of the rf pulse. The line shape of 
the rf spectra has a strong dependence on the rf pulse 
shape. The rf pulse used in Ref. [TtJ is approximately 
rectangular with duration T = 100 /is. We compute the 
experimental rf signal as the convolution of /trap(^) with 
the Fourier spectrum of the rf field [36| 



imation (LDA). In the experiment [l7| not only a single 
(central) 2D layer is populated but also additional pan- 
cakes along the axial direction. For each of these layers 
we calculate the rf response 



^trap 



I 



d'rp^^Ur)ll^\co,r). 



(10) 



trap 



{uj-uj') smc^{uj'T/2). (11) 



The various pancakes are indicated by the index v 
and Plda(^) ~ me^p\r)/27rh'^ is the Thomas- Fermi 
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FIG. 6: (color online). RF spectrum of a trapped quasi-2D 
Fermi gas for (a) B = 227 G, (b) B = 225 G, (c) B = 224.5 
G, and (d) B = 224 G. The attractive polaron gives rise to 
the peak at negative frequencies and the repulsive polaron 
corresponds to the peak at positive frequencies. The dashed 
line depicts the trap averaged signal /trap while the solid line 
shows the expected experimental signal /exp which takes into 



account a rectangular rf pulse of duration T - 
scale) [17.] . 
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In Fig. [6] we show the resulting rf spectra for different val- 
ues of the external magnetic field. Note that the broaden- 
ing may shift the apparent peak position, see e.g. Fig.fTt). 
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FIG. 7: (color online). Same as Fig. [6] but for different rf 
pulse duration (from top to bottom) T = 50 /is, T — 100 /is, 
and T — 200 /xs (not to scale) . 



